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ABSTRACT 



We derive some physical properties of ideal assemblies of identical particles 
obeying generalized exclusion statistics. We discuss fluctuations, and in this con- 
nection point out a fundamental contrast to conventional quantum statistics. We 
demonstrate a duality relating the distribution of particles at statistics g to the 
distribution of holes at statistics 1/g. We suggest applications to Mott insulators. 
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1. Exclusion Statistics 



Several years ago, Haldane [1] introduced the notion of generalized or fractional 
exclusion statistics, interpolating between bosonic and fermionic statistics. Moti- 
vated by the properties of quasiparticles in the fractional quantum Hall effect and 
in one dimensional inverse-square exchange spin chains, he defined the statistics, 
g of a particle by: 



where N is the number of particles and is the dimension of the one-particle 
Hilbert space obtained by holding the coordinates of N — 1 particles fixed. Since 
any number of bosons can occupy a given state, <ijv+AiV — djy, and hence g = 0. 
By contrast, the Pauli exclusion principle implies that g — 1 for fermions. Par- 
ticles with intermediate statistics - Laughlin quasiparticles with g = 1/m and 
1-D spinons with g = 1/2 were the two examples given in [1] - satisfy a general- 
ized exclusion principle. Recently Wu [2] has discussed the statistical mechanics 
of particles obeying a generalized exclusion principle locally in phase space (see 
below). Henceforth we shall call such particles g-ons. Ouvry [3] had previously 
discussed related statistical statistical distributions in the context of anyons in a 
strong magnetic field. Bernard and Wu [4] have shown that excitations in the 
Calogero-Sutherland model obey this statistical mechanics. 

An earlier form of exotic statistics - anyons [5, 6] - has proved very influential 
in the study of two-dimensional systems. Anyons are particles whose wavefunc- 
tions aquire an arbitrary phase e l ° when two of them are braided. Unlike fractional 
exclusion statistics particles, anyons are special to two dimensions since in higher 
dimensions two exchanges can be continuously deformed to no exchange. At first 
glance, exclusion statistics seem to have little to do with the braiding properties of 
particle trajectories which are the starting point for anyons. However in a recent 
important paper, Murthy and Shankar [7] showed that anyons do satisfy a gen- 
eralized exclusion principle (contrary to Haldane [1]). They did this by relating 



9 = 



dN+AN - dN 




AN 
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the exclusion statistics parameter, g, to the high-temperature limit of the second 
virial coefficient, which is non-trivial for an anyon gas [8]. The linchpin of their 
argument is that in a theory with a high-energy cutoff (e.g. any continuum model), 
the transmutation of statistics by attaching flux tubes [6] will generally push some 
states beyond the cutoff, thereby reducing the Hilbert space dimension. This gen- 
erates a fractional exclusion statistics that persists even as the cutoff is taken to 
infinity. 

In this paper, we investigate the elementary statistical mechanics of frac- 
tional exclusion statistics particles further. As an example, we calculate the low- 
temperature specific heat of an ideal gas of ^-ons both for a conserved and a non- 
conserved particle number. In the former case, we find a surprise: the chemical 
potential appears to be independent of temperature to 0(T). We consider fluc- 
tuations in phase space occupancy for g-ons, find that they exhibit subtleties not 
present for bosons or fermions, and on this basis emphasize that one needs much 
stronger assumptions than Haldane's to justify g-on statistical mechanics. We also 
display a duality transformation that relates statistics g particles to statistics 1/g 
holes. Finally, we attempt to make it plausible that electrons in Mott-Hubbard 
systems described by the t — J model might obey g — 2 statistical mechanics; by 
duality, the holes will then obey g — 1/2 statistical mechanics with a modified 
effective temperature. 

We now briefly recall the framework of g-on statistical mechanics [1] [2], to fix 
notations and for later use. The statistical distribution of an ideal gas of g-ons can 
be obtained by maximizing the entropy subject to the constraints of fixed particle 
number and energy, as in [9]. We imagine dividing the one-particle states into a 
large number of cells with I; > 1 states in the each cell and then count the number 
of configurations with ni particles in the i th cell. An elementary combinatorial 
argument gives: 



where d ni is the dimension of the one-particle Hilbert space in the % cell with the 




(1.2) 
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coordinates of ni — 1 of the particles held constant. // one can apply the definition 
(1.1) locally in phase space - a big assumption that we tentatively adopt as the 
working definition of g-ons, but will need to discuss critically and refine shortly - 
then 

d n ,=k- g{m - 1) (1.3) 
Hence we must minimize the quantity, 

S = J2( k + (! - 9)(ni - 1)) ln(A; + (1 - g)rn - 1)) - m Inn* 

i 

- (k - g(rii - 1) - 1) \n(k - g(rii - 1) - 1) - fie^i + Pfirii 



(1.4) 



with respect to the occupation numbers, rif, f3 and are the Lagrange multipliers 
which enforce the constraints of fixed energy and particle number. By standard 
arguments [9] one can show that (3 is the inverse temperature and fi the chemical 
potential. Differentiating with respect to iii, we obtain (neglecting terms of 0(1) 
which are negligible compared to Ink), 

OS 

— = (1 - g) \n(k + (1 - g)rii) -\nm + g \n(k - gm) - (3(ei - /x) (1.5) 
Upon exponentiation this yields, 



rue 



P ^ = (l-gn^(l + (l-g)n^ ° (1.6) 



where we have written n; = ^ for short. This is the fundamental equation that 
implicitly defines the distribution function for g-ons. In general it can only be 
solved numerically, but for special cases including g = 0,1,^,2,^,3,^,4 it can be 
solved analytically. The g = 1/2 distribution function is: 



{ni) g =i/2 = (1-7) 



(l + Ae 2 ^ e ^^y 
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2. Low- Temperature Properties 

At T = the distribution function vanishes for e > /x and takes the value 2 at 
e < pL. More generally, one can show by inspection of (1.6) that at T = 

<*).-{»■ " e< " ; («) 

[0, if e > /i. 

It is quite striking that at T = particles of general exclusion statistics exhibit 
a "Fermi" surface. This fact dictates the low-temperature thermodynamics of 
systems of these particles when the particle number is conserved. 

We now develop an expansion in powers of the temperature - analogous to the 
Sommerfeld expansion for fermions - for the thermodynamic functions of g-ons. 

The Lagrange multipliers (5 and fj, are determined by by the constraints 

N = Y,n l = vJ £^ n g (e(p) ~ u) (2-2) 

/dPp 
(2.3) 

If the single particle energies are of the form e(p) = ap n , then these can be rewrit- 
ten: 

= Jdee^figie-fi) (2.4) 

E f D 

— = dee~ n g (e - fi) (2.5) 

The uninteresting constant factors have been lumped into a. The right hand sides 
of these equations are both of the form 

oo 

J dee p n g (e - \i) (2.6) 
o 

Such integrals are expanded by changing variables to i] = f3e and shifting the 



6 



integral by (3/ji, 



oo oo 

pP +1 JdeePn g (e-fx)= J dr)(r] + ^fn g ft), 



-Pit 

oo 



= J dr)(r) + P/i) p n g ft) + J drift + Pn) p n g ft), 



-fo 
oo 



= J drift + Pn) p n g ft) + J drift + f3^) P {^gft) - ~) + ^ di]-. 

-/fy -/fy 

(2.7) 

Since the integrand of the second integral vanishes exponentially at the lower limit 
of integration, this limit may be taken to — oo at low temperature. This may then 
be combined with the first integral while the third integral may be done, yielding 

oo oo 

(3 p+1 j dee p n g (e - rt) = ± (^) p+1 + j drj^ft + (3f,) p n g ft) + (^ ~ v) P (n g (- V )-^) 



(2.8) 

Expanding in powers of using the binomial theorem, we obtain 

oo 

f3 p+1 J dee p n g (e - = ± (/^ 1 + £( W' Q C i ( 2 - 9 ) 

o i 

where we have isolated the pure numbers 

oo 

Cj = I drirf(n g ft) + (-iy (n g (- V ) - i)) . (2.10) 
o 

Finally, then, we obtain 

,| = ^( 1+9W „) E (I) 3+1 (- 1+ //")q) (2,1) 
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where /iq and E$ are the zero-temperature chemical potential and energy, respec- 
tively. Thus to second order in T 



H = Hq(\ + gCo(-) - g{-\ + D/n)cA— Y) (2.13) 
V/i / V/io^ / 



E — Eq(\+ (<?(! + D/B )0, + g ^( (J/ " >(1 2 +J/n) - (1 + D/nf)) (L 



HO' 

(2.14) 

In the special case g — 1/2, Co can be evaluated analytically. Using the distribution 
function (1.7), one finds rather surprisingly 

oo 

(C ) 0= i/ 2 = / r + / du( r - l) 

9 1 J (l + 4e u )i J Vl + 4e u )3 ' 

-oo v 7 

oo 2 Co 1 K\ 

1 r dv i i r dv , i . 1 ' ; 

~27 V(l + V 2)i + 2J ~\i + v 2)k ~ ' 
2 v ' V 

= 

Hence, the g = 1/2 specific heat to order T is simply 

%) = (l + /n)^l (2.16) 

T/ 5 =l/2 

Ci may be evaluated numerically: Ci = 1.6449. In the case of fermions C n = 
for all even n due to particle-hole symmetry, but this is not the case for g = 
1/2: specifically, Ci = 1.2021. Remarkably, we find numerically that Cq = 
for g = 1/3, 1/4 as well. These results - together with a duality property to be 
demonstrated later which implies that Co = for g = 2,3,4 as well - lead us to 
conjecture that Cq = for arbitrary g. 
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If the particle number is not conserved one has, immediately, 



^-^ e(p)n g (e(p)) = aV T% +1 J d V ^ n 9 ( V ) (2.17) 
For ^-ons with D/n = 1, the integral in (2.17) takes the value 0.9870. 

3. Fluctuations, and a Perspective on the Assumptions 

Let us attempt to find the probabilities for various occupation numbers of a 
single state, for g = 1/2. Defining f(n)e~ n P( e ~^ as the probability for n-fold 
occupancy, we derive from (1.7) the formal relation 

£n/(n)e-^( £ -^ 1 

Zf(n)e- n ^-v) ~ (l +e 2/3( e - M ))i • ( • ' 

Matching coefficients and normalizing /(0) = 1 we find /(l) = 1, /(2) = 1/2, 
/(3) = 1/8, /(4) = 0, /(5) = -1/128, ... . Clearly something has gone awry here! 

The mathematical problem is as follows. The dimension of the many-particle 
Hilbert space 

= (k + { i- g )(N-D)< 

JV!((*-1)- 9 (JV-1))! 1 ' 

for N particles when the cell includes k states vanishes at g — -ffti where m = 
k, k + l,...,k + N — 1. This gives the correct result for fermions, namely W = if 
N > k since g — 1 is then one of the zeroes of W. However, W does not vanish for 
general g when gN > k. Indeed, W can be negative when gN > k. (To see this, 
consider the simple case of g — 1/2 for N > 2k and N an even number.) To get a 
sensible result for the Hilbert space dimension, W, we must stipulate that W = 
if gN > k, which complicates the minimization of the entropy. If we take a large 
cell size k we can safely ignore this complication, since W is small for gN > k. For 
instance, if g = 1/2, N = 2k + 2, W \ and if iV > 2k, then W ^. 
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Thus for g-ons it is important to keep the cell size large, in contrast to the 
Fermi and Bose cases where one could take k — 1 with impunity. In the fermion 
case, for instance, one has 



where pi is the probability of having occupation number i. One can calculate mo- 
ments of the occupation number, n m , by differentiating the distribution function. 
The closure property, n l = n m holds because one can take k = 1 and the only 
allowed occupation numbers are and 1. For ^-ons, one must keep 2k equations 
for n, n 2 , . . . , n 2k in terms of pi,p2, ■ ■ ■ ,P2k- The closure equation which expresses 
n 2k+2 j n terms of the first 2k moments is only satisfied to O(j^). 

One must therefore exercise some care in using the distribution functions de- 
rived from (1.6). The cell size must not become too small in view of the preceding 
paragraph, but it must also not become too large, because if the energy spread 
within one cell becomes comparable to the temperature then the notion of a char- 
acteristic energy for the cell becomes invalid. It is amusing that negative proba- 
bilities appear in this problem in a natural and meaningful way: it is necessary 
that negative probabilities appear in the description of small cells, if independent 
addition of many such is to generate the correct average occupancy for large cells! 

Our mathematical problem reflects a fundamental implicit physical assump- 
tion in the derivation of <?-on statistical mechanics. For bosons or fermions the 
fundamental assumption of symmetry or antisymmetry of the wave-function holds 
rigorously and locally in momentum space. This is enough to allow one to derive 
the appropriate statistical distribution for an ideal gas locally in phase space. To 
derive g-on statistical mechanics as above one must assume that the generalized 
exclusion principle operates on states of nearby energy, and as we have seen one 
must also take a cell size not too small. Without attempting a rigorous discus- 
sion, we can identify qualitatively the physical circumstance under which these 



Po + Pi = 1 



(3.3) 



pi = n = n 2 = . . . = n' 



(3.4) 
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assumptions become plausible. It is that the effective interaction which reduces 
the Hilbert space dimensions should be essentially local in momentum space. Then 
one may apply the counting arguments to cells containing all the states in a small 
range of momenta: this will be a number of states proportional to the volume, all 
with essentially the same energy. 

The known examples of g-ons have this character. Ideal g-on statistical me- 
chanics operates in the Calogero-Sutherland models [4], which feature long-range 
interactions. Anyon models (that is, 2 + 1-dimensional systems with a Chern- 
Simons gauge field) feature interactions which are singular for nearby momenta, 
resulting in a shift of the allowed values of the relative angular momentum between 
two particles, I — >• I + a. Since two particles occupying the same state must have 
vanishing relative angular momentum, an anyon excludes its state from further 
occupation. This is certainly local exclusion, but the second condition is not satis- 
fied: an anyon in any other state also has its relative angular momentum shifted. 
Hence anyons are not ideal g-ons, but interacting g-ons. Whether the ideal g-on 
statistical mechanics provides a useful first approximation in this case, is a question 
needing further investigation. Hard-core bosons on a lattice have g — 1 according 
to Haldane's definition, but are far from being 1-ons according to our definition, 
and their behavior is poorly approximated by fermions: one expects them to Bose 
condense rather than to form a Fermi surface at low temperatures. 

To conclude this discussion, let us finally display the first-order fluctuations 
concretely, for g = 1/2. Using the standard identity relating (An) 2 to the derivative 
of n with respect to the chemical potential, we find: 




(3.5) 



Thus we find that |-ons have sub-Poissonian statistics, as do fermions (An) 2 



n(l — n). In contrast, bosons are super- Poissonian, (An) 2 = n(l + n). 
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4. Duality 



We alluded earlier to a duality property that relates statistics g and 1/g. This 
property may be seen by rewriting the implicit relation for the distribution function 
at statistics g, 

\ 9 l- 
(fH) g ^-^=(l-9(fH) g \ (l + (l-g)(ni) g ) " (4.1) 

in the form, 

(l " 9(fH) g ) = {nifj 9 (l + (1 - .„)(„;)„)' " (4.2) 

or, 

(l - g(ni) 9 ) e-fa-ri = Q - - g (l - 9(fH) g ) ) ^ Q + (l - 1) (l - * * 

(4.3) 

But the implicit relation for the distribution function at statistics 1/g reads 

(*)i/ ff ^ (ei " M) = (l " fag) ^ (l + (l " ^) 1 9 (4-4) 

Hence, 

1 - ^ 5 (/3(e - //)) = -nu g (--(e - //)) (4.5) 
9 9 

This duality relates the distribution of holes in the g-on distribution (where full 
filling is at n g = 1/g) to that of 1/g-ons at g times the temperature, or alternatively 
1/g times the energy and chemical potential. 

This duality is reminiscent of the one found in Chern-Simons models. There 
one describes anyons as charge Q = q objects which acquire a proportional flux 
$ = q/fi, where \i is the Chern-Simons coupling, and thereby have their statistics 
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transmuted by = Q§ = q 2 /^- The fundamental flux tubes then have flux 
$ = 1/q, charge Q = n/q 2 , and the inverse statistics. The thermal duality (4.5) 
is also reminiscent of the g — > 1/g duality in the Calogero- Sutherland models [10]. 
Indeed, for these models thermal duality as discussed here follows from the known 
coupling-constant duality. The concordance of the general thermal duality for 
abstract g-on statistics with the more specific and complete duality for Calogero- 
Sutherland quanta, vividly confirms the idea that these models embody ideal g-on 
statistics. 

5. Remarks on the Mott Problem 

The application of g-on statistics to one dimensional systems, even for the 
soluble models where it is formally correct, is not straightforward. Indeed in one 
space dimension the application of Fermi statistics to derive the low-energy prop- 
erties of systems of fermionic quasiparticles has to be carefully considered, since 
interactions can change the properties qualitatively. The Fermi liquid must be con- 
sidered as one special case of the generic Luttinger liquid [12]. Thus for example 
the thermodynamic properties of edge excitations in the fractional quantized Hall 
states are not correctly reproduced by the ideal g-on formulas, even though the 
electrons are (for example) formally m-ons in the u = 1/m state. Of course the 
bulk filling fraction is nicely consistent with 1/g filling of the magnetic band, but 
this is a much weaker statement. 

More speculative, but if correct probably more useful, is the possibility of 
applications to systems in higher dimensions. For in higher dimensions the phase 
space arguments of Landau apply, as in his justification of Fermi liquid theory, 
and make it plausible that (unlike in one space dimension) the approximation of 
non-interacting quasiparticles is accurate at low temperature. 

There is a class of insulating materials, the Mott insulators, which are anoma- 
lous from the point of view of band or Fermi-liquid theory. They are insulators 
when their valence band is precisely half filled. From the point of view of this 
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paper, it is natural to hypothesize that in these materials the electrons are behav- 
ing as 2-ons [13]. The most important qualitative feature of Mott insulators, that 
is the existence of a gap at exactly half filling, follows directly. Such behavior is 
suggested, but certainly not proved, by the idea (formalized in the t-J model) that 
because of strong on-site repulsion a single electron excludes two states - namely 
states of both spins - from its lattice site. As we have taken pains to emphasize, 
what is needed is local repulsion in momentum space, which could arise directly 
from a long-range force or indirectly through correlation effects. 

In any case, our hypothesis leads to the statistical-mechanical consequences 
derived earlier, which could be tested in experiments or by numerical work on 
models. Most interesting are effects which arise just below half filling. The central 
consequence is the existence at T = of a "Fermi surface" of anomalous size, 
and with anomalous values of the specific heat, Pauli susceptibility, ... . There are 
quantitative anomalies in the experimentally observed normal state specific heat of 
the CuO-based superconductors [14], and in the size of the "Fermi surface" in the t- 
J model as calculated using high-temperature expansions [15] and variationally [16]. 
(In these calculations, the nominal Fermi surface is identified from a strong feature 
in the density-density correlation function.) These anomalies are at least roughly 
consistent with the 2-on hypothesis: the specific heat is substantially larger, and 
the volume of the "Fermi surface" is roughly twice as large, as would be expected 
for ordinary fermions. 
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